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ABSTRACT 

The predictions of general relativity regarding interplanetary time-delay 
measurements are explored in detail. We conclude that a fourth test of 
the theory is now feasible since the modifications of general relativity in¬ 
troduce an extra delay of about 200 psec when radar waves are reflected 
from either Mercury or Venus near superior conjunction. The uses of such 
measurements to investigate the solar corona, a possible time dependence 
of the gravitational constant, and the precession of Mercury's perihelion 
are also discussed. 
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EFFECTS OF GENERAL RELATIVITY 
ON INTERPLANETARY TIME-DELAY MEASUREMENTS 


I. INTRODUCTION 

Although Einstein’s theory of general relativity forms the basis of almost all cosmological 
arguments and has profound philosophical implications, it has been subject to very few experi¬ 
mental tests. The reason is not hard to find. On a laboratory scale, the deviations between the 
Einsteinian and, for example, the Newtonian predictions are almost always too minuscule to be 
detected. In fact, only three tests have been made since Einstein's theory was given its definitive 
form in 1916, all having been suggested in his original papers. Two relate to the interaction of 
matter with electromagnetic waves; the third relates to the interaction of matter with matter. 

The gravitational red-shift experiment, which belongs to the first category, really tests only 
the principle of equivalence and the Einstein prediction is obtained in other theories as well 
(see, for example, Nordstrom ). The most accurate such experiment was performed by Pound 
and Snider who employed the Mossbauer effect to detect the change in frequency of gamma rays 
alternately "rising" and "falling" in the earth's gravitational field. Their results indicate con¬ 
firmation of prediction to about 1 percent. The second test, the prediction of the bending of the 
path of starlight as it passes near the sun, has been subjected to repeated study during solar 
eclipses from 1919 to the present. Such is the difficulty of this experiment, that the results of 
the various observers are consistent with each other and with the predicted value to only about 
25 percent. Most serious reviewers are therefore of the opinion that this prediction of general 
relativity has not yet been verified definitively. The third test involves the prediction that the 
perihelion of Mercury's orbit undergoes precession of 43 seconds of arc per century, in excess 

of the amount calculated from Newtonian theory. The Einstein value seems to have been verified 

4 

to within about 2 percent (see Clemence ). 

A fourth test of general relativity has now been made possible by advances in radar astron¬ 
omy. This test involves measuring the time delays between transmitting radar pulses toward 
either Venus or Mercury and detecting their echoes. These measurements must be taken at dif¬ 
ferent relative orientations of the earth, the sun, and the target planet, with the most crucial 
ones being those near superior conjunction when the radar waves pass closest to the sun. For 
such configurations, as will be shown in Sec. II, predictions based on general relativity indicate 
that the time delays will be increased by as much as 200 psec because of the influence of the sun's 
gravitational field on the speed of radio wave propagation. The increase at inferior conjunction, 
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on the other hand, amounts to only about lOpsec. Hence the difference, which is the significant 
measurable quantity, is almost as large as the maximum value of the increase.* 

The actual test will entail a meticulous comparison of all the observations with the theoretical 
predictions. The unknown parameters (such as the initial conditions of planetary motion, and the 
masses and radii of the planets) will be estimated from the data using the statistical theory of 
parameter estimation.^ The values of the parameters so determined will then be reinserted into 
the theory, and the resulting predictions of time delays will be compared with the observations. 

If the residuals, observed minus theoretical values, are smaller than or comparable to the meas¬ 
urement errors, then we conclude that the experiment "supports" the theory; whereas, if the 
residuals seem to be systematically larger than the estimated errors, then either we have over¬ 
looked some effects on our measurements or the basic theory is inadequate. We could, of course, 
perform the same type of analysis using the Newton instead of the Einstein theory as a basis for 
the comparison. In order to ascertain a priori whether the set of proposed measurements con¬ 
stitutes a test that can distinguish between two theories, we calculate the expected values of the 
measurements on the basis of both, given comparable initial conditions, to determine whether 
the differences exceed the anticipated measurement errors. That the radar experiments being 
proposed here will in fact provide a meaningful test of Einstein's theory can most readily be 
shown analytically by calculating the difference Abetween the proper-time delay predicted in 
general relativity and the corresponding flat-space value. 


II. APPROXIMATE DETERMINATION OF At 

r 

We wish to calculate the proper-time-delay difference Ar^, given the configuration earth- 
sun-planet. First we calculate the delay predicted in general relativity. If the coordinates are 
denoted by x^(p =1—4), then the differential equations for the light rays (which travel along 
geodesic zero lines) are 


,2 p 

d 

, 2 
ds 


+ 


M- 

i P <r 


dx p dx 


a 


g 


ds ds 

dx^ 1 dx 1 ' 
ds ds 


= 0 


= 0 


( 2 . 1 ) 


where the Christoffel symbols 

by 



are given in terms of the metric tensor 


g „(h-. P, <j. V=i -* 4) 


* This relativistic influence on delay had been investigated previously for configurations near inferior conjunction, 
but at that time the delay measurement errors were several orders of magnitude larger than the general relativity 
effect. (In addition, even if the errors were comparable to or somewhat smaller than the change introduced by 
relativity, no experimental test of the theory could have been made. The sensitivity of the time-delay measure¬ 
ments to errors in the Astronomical Unit (AU) and in the radius of Venus cannot be distinguished from the sensi¬ 
tivity of the delays to changes in the relativistic contribution. Even in principle, better than a 1-psec measure¬ 
ment accuracy is required to separate these effects near inferior conjunction. But at that level a whole host of 
obscuring influences come into play, for example, topographical variations on the planetary surface. A realistic 
test of general relativity using only measurements made near inferior conjunction can therefore not be expected 
to be feasible for the foreseeable future.) 
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t We assume, of course, that the data are redundant, i.e., that there are far more measurements than are necessary 
to determine specific values for the parameters. 
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( 2 . 2 ) 


g^ s = g -1 cofactor (g ) 


(2.3) 


and 


P a ’P~ dxP 


(2.4) 


The contravariant metric tensor g^ s is symmetric and satisfies 



a solution to these equations which in empty space is static, has spherical symmetry, and be¬ 
comes the flat metric at infinity. Thus, a spherically symmetric distribution of matter will give 
rise to a gravitational field outside it which is described by Schwarzschild’s solution. We there¬ 
fore use this solution to represent the gravitational field of the sun. 

In order to proceed further, we choose a coordinate system in which to calculate. It might 
be thought that our end result — the proper-time delay between pulse transmission and echo 
reception — would be independent of the choice of coordinate system. However, this conclusion 
is not valid per se because of our assumption that the positions of the earth and planet, relative 
to the sun, are given. Our result will be independent of the coordinate system only if we spec¬ 
ify the planetary positions in an invariant manner. Such a description can be given in terms 
of the line element ds* (The radial coordinate of each body in a given coordinate system can 
be determined from a prescribed value of ds; the angular variables can be treated in a similar 
invariant manner.) But this description would in a sense be begging the question. Our goal is 
really to compare the general relativity prediction with the flat-space prediction, whereas we 
presumably know the planetary positions only in the Newtonian framework. Determining the 
corresponding relativistically invariant positions is not easy. It would probably require a re¬ 
analysis of the optical data upon which the "Newtonian" positions were based. For simplicity, 
and since we are trying only to establish the meaningfulness of the proposed test, we shall choose 
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the usual coordinate system (see, for example, Bergmann ) to represent the Schwarzschild 

solution and shall henceforth ignore the "coordinate-system problem," except for a comment 

7 

preceding Eq.(2.17). In rectangular coordinates, we find 


S44 1 r 


2r 

o 



s 


g = —6 

6 rs rs 


- x x 

r — r A r A s 
o 


r, s = 1 -*-3 


(2.6) 


* This approach was suggested by L. Witten. 
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where 


and 



= c t 

o 


GM 

r = —7T— —1.5 km 
o 2 
c 

o 


(2.7) 


( 2 - 8 ) 


since G is the gravitational constant, M , the mass of the sun, 

s 


and c 

o 


the speed of light. 


3-63 3225 


ACTUAL PATH 



Fig. 1. Geometry of path of radar pulse: 
rectangular coordinates. 


Using Eqs. (2.1) and (2.6), we can find the round-trip coordinate-time delay between pulse 
transmission and echo reception. An approximate result can be obtained almost immediately by 
assuming that the "spatial" path is a straight line (see Fig. 1) and choosing the coordinate axes 
so that this path is parallel to the x-axis. We also assume, for simplicity, that the earth re¬ 
mains stationary between transmission and echo detection. Although affecting the time delay 
slightly, this assumption has a negligible effect on Ar^. (However, see Sec. III.) 

The straight-line path assumption has two effects on the calculated time delay: 

(a) The decrease in path length (see Fig. 1) causes the actual delay to be 
underestimated. 

(b) The closer approach to the sun of the straight-line path results in the 
light ray passing through a higher gravitational potential, and suffering 
therefore a greater decrease in speed and a consequent increase in 
delay. 


We can estimate quantitatively the contribution of each. For the first effect listed, we compare 

the length of the straight-line path with the two straight-line segments, one tangent to the actual 

light-ray path at the earth's position, the other tangent to the actual path at the planet's position, 

as shown in Fig. 1. The mutual inclination of these two segments, which represents the total 

angular deflection of the ray, we call y. As was originally shown by Einstein, even for d equal 

to one solar radius (d = R ), the angle y is no more than 1.75 seconds of arc, i.e., about 8.5 X 10 

s 


-6 


radians. Assuming, again for simplicity, that 
length L' satisfies the inequality 


= x 


it is easy to show that the actual path 


L<L,< ^ “ L ( 1 + ^)” L ( 1 + M ' (29) 

where L is the length (expressed in light time) of the straight path, and e is defined in the fig¬ 
ure. Hence, 
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^ 12 . 

ALiL'-L<L^-<9X10' L 


( 2 . 10 ) 


Even for L = 10 2 sec, we find AL to be only about 9 nsec which is completely negligible consider- 
ing presently achievable accuracies. 

To estimate the contribution of the second effect, we use a crude approximation to the gen¬ 
eral relativistic expression for the speed of a light wave in the presence of a gravitational po¬ 
tential 4>: 


o 


1 + 


2 $ 


o 


( 2 . 11 ) 


where c q is its speed in the absence of a gravitational field. The difference Ac between the speed 
of a light wave at the point A of closest approach to the sun along the straight-line path and the 
speed at the intersection A' of the two line segments that are tangent to the actual path at earth 
and planet, respectively, is then given by 

Ac = c(A) - c(A') = ~ (4>(A) - <f>(A 1 )] 

o 


2GM s 


c 

o 



1 

(d Q + Ad) 


2GM 


s Ad 


c 

o 



( 2 . 12 ) 


where 


Ad 


L 

4 


sin e 



(2.13) 


Equation (2.12) represents an upper bound on the maximum difference in speed along the two 
paths. Hence, a gross upper bound on the effect AL on time delay of assuming a straight-line 
path is given by 


AL < L 


Ac 


c 

o 




(2.14) 


For L = 10 3 sec, we find AL < 2 X 10" sec even when d = R ^ « 2.3 sec. Hence, this effect of 

o s 

the approximation is also negligible. 

This type of result is actually a general one for refraction phenomena where the change in 
index is small: The effect of the change in path on the time delay is of higher order in the change 
of index than is the effect of the change in speed. 

Having justified our straight-line path assumption, we now calculate the coordinate-time 
delay t explicitly. Since we have determined the path, albeit by assumption, we need only the 
second of Eqs. (2.1) to determine t. Considering the rectilinear nature of the path and its direc¬ 
tion parallel to the x-axis, we find from Eqs. (2.6) and (2.7) that 


2,2 , 

g 44 C o dt +S 
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dx 2 = 0 


i.e., that 
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dx 

c " dt “ c 




/_ °44 \ 

o\ g H / 


= c 


o 


2r 


1 - 


o 


1 + 


2r 2 
o x 


r — r 2 

or- 1 


1/2 


o 


4?(-4)l 


(2.15) 


Since t increases monotonically over the round trip, whereas x first decreases monotonically 
and then increases monotonically, we find that the round-trip delay t is given by 


t = -fp^ + r 

r Jx c J x 


e dx 


= 2 


.S' e <r J e 1 + ~f (* + h)\ dx • (2 - 16) 


P 


P 
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Using Dwight leads to 


? 2r 
t = — (x — x ) + —- 
r c e n c 
o o 


( x + r v /X x\i 

X n + r n)"( 1 V ” ^)] ' 


x + r 
P P 


(2.17) 


where in consideration of our straight-line path. 


2 2.2 

r = x + d 
e,p e,p o 


(2.18) 


What we seek, of course, is not the coordinate-time delay, but the proper-time delay t 
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which is measured by the earth observer. As is well known. 


r = — p ds 
r c 

o ^s. 

1 


(2.19) 


where s^ is the (four-dimensional) position of the earth at transmission, and s ? is the corre¬ 
sponding position at echo reception. Since, for the earth, the geodesic is given by 


ds = g dx^dx^ , 


( 2 . 20 ) 


and since we assumed that the earth remains stationary between pulse transmission and echo 
reception, we find 


ds2 = g 44 c o dt2 


( 2 . 21 ) 


and hence,^ 


f l 2 1 

' 1 g 4 


^dt» « 2 -1 


•'(*-£) 


( 2 . 22 ) 


But (t^ — t j) is just the left side of Eq. (2.16), and we obtain, finally. 
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* If a harmonic coordinate system were employed (see Fock ), Eq.(2. 17) would not contain the nonlogarithmic 
term inside the brace. ^ However, in the harmonic coordinate system the numerical values for the planetary 
positions are related to the relevant invariant quantities in a different way than in our coordinate system. As we 
explained earlier, these differences are of no intrinsic importance. 

t The contribution of the earth's orbital motion to the ratio ~ *l) ' s actually r Q /(2r e ). But this contribution 

is the same as the corresponding contribution to the flat-space delay. Since we are here concerned only with 
the difference in these delays, neglect of the earth's motion is not serious. 
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2r 

T « - (X — X ) + -- 

r C 6 D c 
o ^ o 


2 log 


x , 2 x , 2 . 1/2 
x + (x + d ) ' 
e e o 

e , , 2 ^ , 2,1/2 
x + x + d ) 
p e o 


2x_ — x 

2 p 


X 


L(x 2 + d 2 ) 1/2 
L e o 


, 2 , . 2 , 1/2 
x + d ) 
p o 


(2.23) 


The calculation of the corresponding flat-space delay t, is trivial; the result is simply 

2 


To = - (X — X ) 

fs c e p 
o 


(2.24) 


and hence the difference A t^ is given by 


2r 


2 ± , 2.1/2 
x + (x + d ) 
o u , e e o' 

AT r*— 21o *e- + , 2 + .2.1/2 

o x + (x + d ) 

P P o 


2x — x 
e P 

, 2 x , 2 , 1/2 
(x + d ) 
e o 


(x 2 + d 2 ) 1/2 
p o 


(2.25) 


This expression simplifies considerably in several cases. Near superior conjunction, for 
example, we find* 


4r 

• 

4x x 

y 3x — x v 

A o 

At ~ - 

r c 

o 

io g e 

e p 
d 2 

- ( K p ) 


o 


and near inferior conjunction. 



4r 


X V 

/X — X V 1 

At « — — 
r c 

lo ^e ( 

y H 

* 

o 


k P 7 

\ e ' J 

At elongation, Eq.(2.25) reduces to 


4r 

1 

2x , 

1 

A o 

At «- 

r c 

o 

lo ®e ( 

,ir - 

o ' 

1 ; x « 

J p 


d « x , x 
o e 1 p 


d « x , x 
o e p 


, 2 2 
d « x 
o e 


( 2 . 26 ) 


(2.27) 


(2.28) 


This last form is valid only for Mercury since, for Venus, x g » d Q at elongation. 

Although g and g contribute equally to the dominant logarithmic term in the expression 
44 r b 

for At , their relative contributions to the second term vary with the path of the radar pulse, 
r 

The maximum magnitude of the difference between these contributions is, however, too small 
to be reliably detected experimentally in the near future. 

In Fig. 2, we show the values of At^ as a function of the angle between the earth-sun line 
and the earth-Mercury line. This function, of course, is double-valued, one branch corre¬ 
sponding to the planet's being on the same "side" of the sun as the earth, the other branch to the 
planet's being on the opposite side. Both branches are shown. Figure 3 presents the same re¬ 
sults for Venus. In constructing these figures the orbital eccentricities were neglected, and the 
round-trip delays corresponding to several of the orientations were included. 


III. RIGOROUS DETERMINATION OF A 

We now proceed to a rigorous derivation of Eq. (2.25). That is, here we shall make a formal 
expansion of the solution to the equations of motion in powers of r . Using Eqs. (2.1) and (2.6), 


2 2 1/2 2 21/2 2 
Note the identity: -log [-x + (x + d ) / ] = log [x + (x +d) 1 - log d 

e e e 


7 































A r r (psec) 





Fig. 2. Effect of general relativity 
on earth-Mercury time delays. 


Fig. 3. Effect of general relativity 
on earth-Venus time delays. 
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Fig. 4. Geometry of path of radar pulse: 
polar coordinates. 
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TWO-WAY DELAY (AU) 
























11 

and making the transformation from rectangular to polar coordinates, we find the equations 
of motion of a light ray to be 


< a |) 2 - e '^) 2 - r 2 (^) 2 = 0 


= k 


ds 


di 


ds 


2d£ = 
ds 


h 


(3.1) 


where r and cp are polar coordinates in the plane of motion of the light ray, and t is the time 
coordinate. The quantity e^ stands for 


2r 

u A ° 

e r = 1- 


(3.2) 


whereas k and h are constants of integration with s being a parameter. 

Our goal can be accomplished as follows: We assume as before that the positions of the earth 
and target planet are given and that the earth remains stationary between transmission and detec¬ 
tion of the radar pulse, and we calculate the round-trip time delay twice, once using the values 
for g v given in Eq. (2.6), and once letting g = ±6 with the plus sign holding for the temporal 
component and the minus sign for the spatial components. Each of these values for the difference 
in coordinate time is then separately converted to the appropriate proper-time delay. The two 
values are then subtracted, yielding At . One aspect of the subtraction operation deserves 
special mention: The distance of closest approach is different in the two calculations and the 
relation between them must be established to the proper accuracy. 

Let us do the flat-space calculation first. We refer to Fig. 4 for the appropriate geometry. 
For convenience, we orient the x-y axes such that, for the flat-space path, the point on the path 
closest to the sun lies on the y-axis. (This orientation, as before, introduces certain algebraic 
simplifications in the calculations.) We denote this distance of closest approach by d Q in the flat- 
space case, and by d in the curved-space calculation. 

For the flat-space case, of course, we could deduce the result by inspection. However, to 
illustrate the general method we proceed as follows: Here Eqs. (3.1) reduce to 

. dt .2 , dr .2 2. d (p .2 n 

( di' -- r = 0 . 


ds 


ds 


dt 

ds 


= k 


2 dip 


ds 


= h 


(3.3) 


Combining the second and third yields 


dt _ 2 k. 

d (p 1 h 


(3.4) 


2 

whereas multiplying the first by [ds/d(p) leads to 




(3.5) 


Throughout we use units in which is unity. 
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Substituting (3.4) into (3.5) yields 


i.e.. 


4k , dr,2 2 n 

r - r *° 


. dr ,2 k 4 2 

*5?' "72 r “ r 


(3.6) 


Defining 


R _ h 

R = k 


(3.7) 


we obtain 


d cp = ± 


dr 


<?-*') 


1/2 


= ± 


Rdr 


, 2 D 2 , 1/2 

r(r — R ) 


(3.8) 


Since cp is monotonically increasing for the first half of the round trip and monotonically de¬ 
creasing for the second half, whereas r first decreases to d Q , then increases to r , etc., we 
see that the minus sign applies for the first part of the one-way trip and the plus sign for the 
second, etc. We will calculate explicitly only the one-way delay: The two paths are identical 

o 

so that the total delay t^ g is just twice the one-way result. Integrating (3.8) yields 


, -1 , R v 

<p = ± COS ( — ) + <p. 


n ^ -1 . R . . 7T 

0 ^ cos (7)^2 


(3.9) 


where the minus sign holds for 0 (p ^ 7r/2 and the plus sign for 7r/2 <</?<: 7r 
R and cp ^ in terms of the boundary conditions: 


We can determine 


cp - —cos 


cp = cos 
P 




+ <P 1 


These equations can be solved in a straightforward manner to obtain R and cp^. 
the first from the second and taking the cosine of both sides, we obtain 

/ x 1 rD 2 , 2 0 2,l/2 , 2 0 2. l/2, 

cos ( (p — cp ) - - [R — (r — R ) ' (r — R ) 7 1 

p ^e r r 1 e p 1 

^ e p ^ 


(3.10) 
By subtracting 


Eliminating the radicals, we find 


[r r cos(v - cp ) — R 2 ] 2 = (r 2 — R 2 ) (r 2 — R 2 ) 
1 e p ^p ^e 1 v e P 


Hence, 


R = 


r r sin {cp — cp ) 

0 p' r p ■ e 


r 2 , 2 0 , ,.1/2 

fr + r — 2r r cos {cp — cp ) 

1 e p e p p ^e ,J 


(3.11) 
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The value of ip ^ then follows directly from Eqs. (3.10). The denominator in (3.11) is simply 


r cos (p + r cos (7 t — ip ) = r cos ip — r cos ip 

0 T 0 p ' • p f 0 T 0 p r p 


whereas the numerator is 


since 


Hence, 


d (r cos (p — r cos cp ) , 

o e e p ^p 


r sin (p = r sin ip - d 
e e p ^ p o 


R = d 


o 


(3.12) 


and inserting this result in Eq. (3.10), we find 


7T 

= 2 


(3.13) 


“1 ”1 — 1 “ 1 

since sin x + cos x = 7r/2, when x < 7r/2; and sin x — cos x = 71 -/ 2 , when tt/2 < x 7 r. These 

results also follow from the fact that 


dr 

dip 


= 0 


r=d 


o 


Thus, 


dr 

dip 


r=d 


d (d 2 - R 2 ) 1 / 2 
o o___ 

R 


= 0 


> R = d 


o 


o 


Since <p - n/z when r = d , by our choice of axes, we have 


7T 


-1 


7 r 


2 cos ( 1 ) + ^ 2 


To obtain the time delay t^ explicitly, we use Eqs. (3.4) and (3.5): 


,dr x 2 ,dip K Z , dr ,2 R^ ( Z\ A 

<dt> = <dT > <d^> = ^4 (^2 ~ r ) = l ~ 


R‘ 


(3.14) 


Using Eq. (3.12), we find 

dt = ± — 


rdr 

. 2 .2.1/2 

o 


(3.15) 


Again, we find t increases monotonically, whereas r decreases monotonically from r e to d^ 
and then increases monotonically to r . The plus sign therefore holds for the second part and 
the minus sign for the first. Remembering the factor of two needed to convert to the round- 
trip value, we find 8 
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'fs 


* 2 i. 


rdr 


, ,2 , 2 . 1/2 
d (r — d ) 
o o 


+ 2 


T p rdr 

K ,2 ,2 

^d (r — d 
o o 


d 2 ) 1 / 2 

o 


. r , 2 , 2 . 1/2 , 2 , 2 . 1 / 2 . ... 

= 2 [(r e _d o ) -< r p - d 0 ) ) > < 3 - 16 ) 


But, as mentioned above, 


d = r sin (p = r sin (p , 
o e ^e p p 

and hence, reinserting explicitly, we find 


(3.17) 


t r = — (x — x ) 
fs c e p 
o 


(3.18) 


since x is negative for the configuration of Fig. 4. 

I 

Now we proceed analogously for the curved-space case 
members of Eqs. (3.1), we find the analog to Eq. (3.4): 

dt -Li 2 k 
d (p h 


Combining the second and third 


(3.19) 


Multiplying the first by (ds/d (p) , as before, leads to 


u.< dt x 2 -Li. dr ,2 2 n 

e W " e ( d^> " r = ° 


Substituting (3.19) into (3.20) yields 


(3.20) 


i.e.. 


"P r 4 iL_ _ e "P/_dr \ 2 _ r 2 - 0 

r 2 e 1 dcp ' r 




4 k 2 u 
= r —y — r e r 

h 


(3.21) 


Again defining R as in Eq. (3.7), we obtain 

, , _ Rdr _ 

d(p ,2 p 2 p.l/2 

r(r — R e r ) ' 


(3.22) 


As before, we note that for the configuration drawn in Fig. 4, <p increases as r first decreases 

from r to d and then increases from d to r . Then (p decreases as r first decreases from 
e p 

r to d and then increases from d to r . Since the round-trip is symmetric, we will consider 
P e 

only the one-way case and multiply the result by two. We have, therefore. 


<p = - 



Rdr 

-rV ) i/2 


r g ^r^d ; (P Q 4(p4<P( d) 


Rdr 


r(r - R eH 


2 m 1/2 


d ^ r ^ r p ; </?(d) 4 (p 4 <P p 


(3.23) 


Let us first consider the case r g ^r ^d. By changing the integration from r decreasing to r 
increasing, we have 
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But 


<P = 


f 



Rdr _ 

- R 2 e M ’) 1//2 


2 r 

u . o 
e n = 1- 


and since r = d at (dr/d (p) - 0, we find from Eq. (3.21) that 

,2 / 2r 


R‘ 


-(*-? - • 


i.e., that 


a 2 , -i 


2 r 


1 - 


o 


or 


R 


(*- 


_o \ v o' 

V 


(3.24) 


We determine d from the initial conditions, after integrating Eq. (3.23). For this latter purpose, 
we substitute Eqs.(3.2) and (3.24) into (3.23): 

dr 


<P = 


/ 2r \ 

(“T 2 


•/2 J 


2 

r — 


/ 2r \ 
*-T 


/ 217 \ 

-- 


1/2 


= J „l/2 r „3. .3 


dr 


r A/ ^ [r"d - d"r - 2r (r 3 - d 3 )] 1 ^ 2 

l Q j 


(3.25) 


But 


Hence, 


r 3 — d 3 = (r — d) (r 2 + rd + d 2 ) 


- j3/2 i 


cp = d 


dr 


[r(r — d)] 1 / 2 [rd(r + d) — 2r^(r 2 + rd + d 2 )] 1 ^^ 


(3.26) 


Since r Q « d, we expand: 


dr 


^ ~ ^ \ z 2l/2 
J r(r — 6T) /L 


r r d 

1 + -£ + 0 


d r(r + d) 


and also find 


8 


(3.27) 


dr 1 -1 IdI 

---~ = -t cos 

. 2 ,2.1/2 d 'r 1 

r r — d ) ' 


A . -1 I d I . 7T 

0 ^ COS |-| 4 2 


(3.28) 
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To evaluate the second integral, we note that 


1 


r 2 (r + d) dr 2 


+ 


d 2 r d 2 (r + d) 


(3.29) 


Hence, 


<P 


(l + -f ) cos' 1 |f | + r o d 2 j 


dr 

Tl , 2 , 1/2 

(r - d ) 7 


1 


+ 


dr 2 d 2 r d 2 (r + d) 


8 ^2 

Using Dwight and Grobner and Hofreiter, we find 


. (3.30) 


Cp = cos 


-1 i d | , o,2 j2 x l/2 r 2r + d , , 
1 r 1 d l r(r + d ) 1 


(3.31) 


We could determine the constants of integration, d and (p , , in terms of the boundary conditions. 
Equivalently, let cp d be the angle at which r = d. Then, we have 


<? d ~ = COS 


“ <P d = cos 


d 

^ r o , 2 ,2.1/2 

+ ^ (r e “ d > 

r 2r + d 
e 

r 

e 

r (r + d) 
L e e 

d 

+ fo (r 2 _ d 2 )l/ 2 

d ' p ' 

r 2r + d 

P 

r 

P 

r (r + d) 
1 P P 


(3.32) 


Since we are solving only to first order in r Q , we can replace d by d Q in the coefficients of r^. 
Adding the two members of Eq. (3.32), we find 


cos 


/ -1 I d 1 , -lid.\ 1 r ,2 .2 ,2,1/2. 2 ,2,1/2, 

(cos — + cos — ) = - fd — r — d ) ' r — d ) ‘ 1 

\ 1 r 1 1 r 1 / rr L ' p o ' e o 1 


e p 


= cos ((p — cp ) + f sin [(p — cp ) , 

XY p e d Q VY p v e 


(3.33) 


where 


(r 2 — d 2 ) 1//2 (2r + d ) (r 2 - d 2 ) (2r +d) 

f _ p o _ p o e o _ e o 

r (r + dj r Jr + dJ 


P P o 

From Eq. (3.33) it follows that 


e e o 


d = 


r r {1 — 
e p 


o 


cos ( cp — cp ) + — 7 — f sin {ip — <p ) 
^p e d^ ^p e 


2 | 1/2 


2 x 2 •> 

r + r — 2r r 
e p e p 


cos (cp —</?)+ 

^ d Y e d 


o 


o 


f sin (cp — <p ) 
^p e 


1/2 


(3.34) 


(3.35) 


which is the analog of Eq. (3.11). Expanding the right side of (3.35) leads to 


r r sin (c p — cp ) 
e p ^ p ^ e 


r 2 i 2 , . ,, 1/2 

t r e + r p - 2 r e r p cos(^. p -^ e )] 


1 - 


r f 
o 

d 

o 


ctn (cp — (p ) -~~ 


r r sin (cp — cp ) 

0 p_ y p e 


r + r — 2r r cos {cp — cp ) 
e p 0 p ~ p t e 


(3.36) 
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Using Eqs. (3.11) and (3.12), we find 

r 

o 


d « d I 1 + -r~ f 
o\ d 

o 


o 


r 2 , 2 - . .. 1/2 

fr + r — 2r r cos (<p — ip )] ' 

1 e p e p xv p J 


ctn ((p — <p ) 

^ p ^ e 


But from the law of cosines, we have 


r 2 , 2 , , .,1/2 

[r +r — 2r r cos (<p — <p )] ' = x — x 
1 e p e p ^p ^e J e p 


It also follows from the above that 


and that 


Hence, 


d + x x 
o e p 


cos (<p — <p ) = 

p ' e r r 


sin (<p — m ) = 

p Y e 


e p 


d (x — x ) 
o e p 

r r 
e p 


(3.37) 


(3.38) 


o 


X — X 

e p 


d + x x 
o e p 

r r 
e p 


r r 
e P 


d (x — x ) 
o' e p 


x x 


(x -x ) 
e p 


which implies that 

d ~ d 


o 


o 



r 

, X X 

V 1 

1 - 

o 

/ e p 

\ f 

, 2 

lx — X 



d 

o 

v e p 

/ J 

and 

X 

< 0. But 

we 


P 




(3.39) 


is not obviously in error. To find we can use either of Eqs. (3.32). Let us take the first: 

d 

< 

r 


(p d = + cos 


-1 


g--L- ( Vp \ f l + ^o (r Z_ d 2 1/2 

dr Vx — x / d v e o 

e o e v e p' J o 


2 r + d 
e o 

r (r + d ) 
e v e o 


We note that if 


-1 A 

cos A = B 


and 


-1 


cos (A + e) = B + 6 
then for e, 6 « 1, 

A + e = cos (B + 6) « cos B — 6 sin B , 

and 

5 - cos B — (A + c) _ c 

sinB ~ sinB 


In our case. 


B = cos 


*&) 
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Hence, 


d 

cos B = —— 
r 

e 


sin B = — 
r 

e 


and 


cos 


l e oe'e p' J x e' o < 


But, from Eqs. (3.10), (3.12), and (3.13), we find 

d 


cos 




<p. 


xx 1 ,r 

e - P . I f 


(x — X ) 
e p 


k) 


and therefore, 


<P t 


/ x f V 

r 

( p ) 

+ x 

V X — x / 

d e 

v e p' 

o 


2r + d 
e o 

r (r + d ) 
e e o 


Using the expression for f from Eq. (3.34), we get 

f = - 


x (2r + d ) x (2r + d ) 

p p o . e e o 


r(r + d ) + r(r + d ) 

p p o e e o 


and 


v o . X -1 

*d = 2 - d” (x e - x p> 
o 


x ^(2r + d ) x ^(2r + d ) 

p p o e e o 


r (r + d ) 
P P o 


r (r + d ) 
e e o 


(3.40) 


Is this result reasonable? If r = r , then x = —x by our original assumption, and therefore 

^ ^ P 

ip^ = 7r/2 as we would expect from symmetry arguments. If |x | « x^, we find <p^ > tt/ 2. This 
deduction is also in accordance with expectations (or so it seems on cursory inspection!). At 
least (3.40) has the expected symmetry properties under reflection through the y-axis. 

To obtain the time delay, we proceed as before. From Eqs. (3 19) and (3.21) we find 

d(? 2 , dr\2 dr 2 2p Fti / r^ _ 2 p\ 

( ~3t ' ( dtp ’ - ( dt' “ 6 „4 \ r 2 r e ) ’ 


and 


dr p/, R 2 p\ 1//2 

ns - 7 ’ ) 


(3.41) 


Thus, the coordinate delay from r = r g to r = d, is 


t 


ed 



and the corresponding delay t 


dp 


from r = d to r = 


r is 
P 



(3.42) 


(3.43) 
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To evaluate the indefinite integral to first order in r^, we find, using Eqs. (3.2) and (3.24), that 


(* - ** •-) 


TJz 


. 2 . 2 .- 1/2 

(r - d ) 




Hence, 


ed 


■I 


e rdr 

, 2 ,2.1/2 r o 

d (r - d ) ' 


c 


2dr 


2 2 
*- ; d (r — d ) 
L o x o 


!/2 + I 


d dr 
Q 


, , x . , , 2 . 2 . 1/2 

d (r + d ) (r — d ) ' 
o o o 


(3.44) 


Integrating yields 


8,12 


. , 2 , 2 . 1/2 , 

ed e o 


, , 2 , 2 . 1/2 , 2 , 2 . 1 / 2 , 

r + (r — d ) ' (r — d ) ' 

-> i e e o . e o 

2 log e d_ (r + d ) 


o 


e o 


Thus, the total two-way delay t is given by 


t E 2 t , + t , = - 

r ed dp c 
^ o 


. 2 , 2 . 1/2 , . 2 , 2 . 1/2 

(r e - d ) + (r - d ) 


4r 


o 


o 


iog. 


r + (r — d ) 
e ' e o 


2.1/2 


o 


r + (r — d ) 
P P o' 

-a 


2.1/2 


o 




. 2 . 2 , 1/2 , 2 , 2 , 1/2 
(r — d ) 1 (r — d ) ' 
e o p o 

(r + d ) (r_ + d_) 


e o 


P o 


(3.45) 


where we have reinserted explicitly. In view of Eq. (3.39), we can write: 



which leads to 


Similarly, 


. 2 . 2 , 1/2 
(r — d ) ' 
e 


1 + r 


o x (x — x ) 
e e p 



- d 2 ) 1 / 2 



1 + r 


o x (x — x ) 

P e p 


We can therefore rewrite Eq. (3.45) as 


(3.46) 


(3.47) 


4r 


t = — (x — x ) + 
r c e d 
o 


o 


o 


( r+Xv/r— x\ . ✓ x Xvi 


(3.48) 


which is exact up to terms of order ( r 0 / c 0 ) • 

Both Eqs. (3.18) and (3.45) represent the time delays in coordinate time. However, as be¬ 
fore, we seek the proper time measured by the earth observer. We therefore use Eq. (2.19) to 
convert each of these results to proper time. Consider Eq. (3.18) first. Here Eq. (2.20) becomes 
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(3.49) 


2 2 2 
ds = c dt 
o 

2 2 

since we have assumed that the earth remains fixed which implies that (dr) = (dc/>) 
the nonflat metric, we have 


= 0. For 


ds2 = °o (’ - Tr) 


dt 


since the temporal part of the metric depends on the gravitational field. Thus 

2 r »l /2 


nt / 2 r . 1/2 , r . 

■ 'M ’ 

«^o ' e ' x e ' 


(3.50) 


(3.51) 


and to first order in r , we find 

o 


4r 


2 . \ , o 

t = — (x — x ) + - 

r c e d c 
o v o 


( ? r + x v /r — Xv ,/X x x — x 

_e_ e\ { P P\ 1 / e _ p x e p 

d ; \—a —/ 2\— — + ~r 

o ' ' o 7 v e p e 


)l ■ 


whereas 


(3.52) 


T~ = - (X — X ) 

fs c o e p' 


(3.53) 


Since 


l0 §( i -r i ) = l0 g ( FT^ ) * 


we can write 


2 r 


At = t — t~ « — 
r r fs c 


o 


o 


( r +Xv / ^x — x Xvi 

r n +X n) _ ( *7^ _T ^)] 


2 x — x x 
e p 


r + 

P P 


(3.54) 


The above formula was derived on the assumption that the planet and earth are on opposite 
sides of the y-axis. Suppose both are on the same side; then the left side of the second equation 
of (3.32) should read (p ^ — <p^. Hence 


cos 


/ -1 I d , -1 I d A 1 . ,2 , , 2 ,2.1/2 , 2 .2.1/2. 

(cos |—| - cos I — |j = —- [d + (r -d ) (r - d ) ] 

' e p' ep ^ 


= cos (ip — cp ) + - 3 — f* sin (<p — <p ) , 

^ p ^ e d ^ p ^ e 

1 o ^ 


(3.55) 


where 


f* = 


, 2 , 2 .l/ 2 .- , x 

(r — d ) ' ( 2 r + d ) 

e o e o 


r (r + d ) 
e e o 


.2 , 2 . 1 / 2 ,, , . . 

(r — d ) ( 2 r + d ) 

p o p o 

r (r + d ) 

P P o 


As previously, we find 

d ’ « d 


o 


r 

/XX, 1 

1 --5- 

e P \ f. 

, 2 

lx - X ) 1 

d 

A 

' e p/ J 


(3.56) 


(3.57) 
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We also obtain the analog of Eq. (3.45) as: 


V = 2(t , - t , ) = 
r ed dp c 


2 


,, 2 ,2,1/2 , 2 ,2,l/2i . 4r o 

I (r — d ) 7 — (r — d ) 7 


o 


P 


1 4 


o 


r +(r 2 -d 2 ) 1//2 

^S e ~ -1-2^72 + ! 

r + (r -d ) ' 

P p o 


/ 2 , 2,l/2 

(r — d ) 7 
e o 

(r + d ) 
e o 


, 2 aVAn 

(r p ~ d Q ) 

(r_ + d ) 


o 


(3.58) 


But, in analogy with Eqs. (3.46) and (3.47), we find 


, 2 , 2 , 1/2 , 2 , 2 , 1/2 

( r e - d ) - (r - d ) ' 


Hence, 


4r 


t' = — (x — x ) 4 
r c e p 
o 


o 


o 


~ x — x — r f' 
e p o 


z r + x v . / x x v i 
log e(r + x )- z(t~ ~ r?)] 


r 4 x 
P P 


and finally. 


4r 

r / r 4 x x 


z 2x — x 

X V 1 

A . o 

At 1 = - 

iog (—■ -■ 

- | 

e 

p - p ) 

r c 

o 

e v r p + y 

2 

r e 

P / J 


(3.59) 


Comparing (3.59) and (3.54), we find that, in general, 

4r 


At = 
r 


o 


o 


. r +Xv . / ^x — x Xvi 

E -r E ) 

' n o' ' e n' J 


2x — x 
e p 


r + x 
P P 


+ O 


o 


(3.60) 


Considering Eq. (2.18), we see that the result (3.60) is in precise accord with Eq. (2.25). 

At this point we should mention that in an interplanetary radar experiment, in addition to the 
time delay, the Doppler shift of the radar wave is also measured; but, although the effect on 
time delay of the change in c is cumulative, the corresponding general relativistic effect on 
Doppler cancels out over the round trip.' In our development thus far we have assumed that the 
earth remains stationary. Whereas such an approximation has no significant effect on At , the 
change Af^ in Doppler attributable to At comes about solely because of the time variation of 
At^ introduced by the relative motions of earth and planet. Since, to first order in v/c, the 
Doppler shift Af is related to the transmitted frequency f and to the time delay t by 

Af = -fT , (3.61) 

and since At^ is most sensitive to changes in d , we can approximate Af^ by 


Af 


/3At \ 

/ 9d \ 


1 r 1 

I °l 


\9d o / 

\ at / 

} 


(3.62) 


where by straightforward differentiation of Eqs. (3.60) and (3.11), we find 


9At 


4r 


o 


8d 


o 


o 


o 


_i_S_ , « pe ~ X P X P 

r (r + x ) r (r + x ) 2\ 3 3 

0' 0 0' P'P p 7 \ r r 


(3.63) 


* If a suitably calibrated frequency-measuring device were stationed on one of the inner planets or in orbit in 
their vicinity, then the general relativistic effect on frequency, the so-called gravitational red shift, would 
indeed be detectable, amounting to a change in frequency at X-band of about 20 cps at the orbit of Venus. Of 
course, when the receiver is inside earth's orbit, the change is a violet shift. 
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and 


9d 


o 


Q . = r r ((D — (p ) 

at e p p e 


cos ((D — (D ) 

^ p r e 


r 2 , 2 0 , \A/2 

f r e +r p -Z r e r p co B^p-V e )] 


• 2 / X 

r r sin ((O — cp ) 
e p e 


r 2 , 2 0 N ,3/2 

t r e +r p - 2r e r p cos(^ p -<p e )] 


Near superior conjunction, this approximation to Af^ reduces to 


(3.64) 


8r 

Af * ± f —- 
r c 

o 


x e x p (<P 

d (x 
o e 



) 


(3.65) 


where the minus sign is to be used for the pre-superior-conjunction configuration and the plus 
sign for the post-conjunction configuration. Even at X-band frequencies and with grazing in¬ 
cidence (d ~ R J, Eq. (3.65) yields only a 3.6-cps effect. Furthermore, since the planetary 

O b 

orbits are not coplanar, the magnitude of ad ^/at will in general be less than indicated by 

Eq. (3.64). At d = 3R , the closest approach possible with the Haystack antenna (see Sec. IV), 

o s 

Af would fall to about 1 cps or less. Aside from the questions of frequency stability and of 
r 5 

other influences on the Doppler shift, it is doubtful that the center frequency of the Doppler- 
broadened echo from a rotating planet could be located with sufficient accuracy to detect Af^ 
reliably. 


IV. OTHER INFLUENCES ON TIME DELAYS 

Are the relativistic effects on interplanetary time delays likely to be obscured by others? 

The most important candidates in this latter category are the imprecise knowledge of plan¬ 
etary orbits and the presence of interplanetary plasma. A moment's reflection suffices to 
show that the orbits of the earth and the target planet can be determined with more than the 
required precision from optical and time-delay measurements distributed around the orbits of 
both planets. For example, time-delay observations of Mercury from earth could be made at 
all positions of Mercury along its orbit with the radar wave never passing near the sun. Ob¬ 
servations along one half could be obtained from elongation to elongation through inferior conjunc¬ 
tion, and along the other half during a corresponding elongation-inferior conjunction-elongation 
period, but with the earth on the opposite side of its orbit. In addition, since the received power 
varies with the inverse fourth power of the interplanetary range, measurements made near in¬ 
ferior conjunction will undoubtedly be more accurate than the corresponding ones made near 
superior conjunction. Hence, the precision of orbital determination should be at least as high 
as the accuracy with which one can make the crucial time-delay measurements near superior 
conjunction. More generally, we can see from this analysis that the sensitivity of the time delays 
to changes in Awill be different from the corresponding sensitivities to changes in the initial 
conditions of the orbits and in the planetary masses and radii. A parameter characterizing 
At^ could therefore be estimated from the data simultaneously with the other relevant ones, 
without incurring any appreciable accuracy penalty from inseparability of effects. The topo¬ 
graphical variations on the target planets are probably small enough so that even the most ac¬ 
curate measurements will not be significantly degraded thereby. 
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The effect At^ of the interplanetary medium on the two-way time delay can be represented 


by 


At « - -■ = X - 1 -— f e N(f) d t sec 
m ) 

f c 

o p 


(4.1) 


where N is expressed in electron/cm , f in cps, c in cm/sec, and £ in cm. Equation (4.1) is 

valid only if f lies sufficiently above the plasma frequency of the medium; this condition will 

hold for all our considerations. Using recently compiled results on the solar corona (see Fig. 11 ;e 
13 

in Erickson ), we find that during a "quiet-sun” period, N(r) can be represented reasonably 
well by 


N(r) =5X10 


(5?) 2 «/cm : 



(4.2) 


from r = 4R to r = 20R . Inside this range the actual N increases more rapidly with decreasing 
s s 

r, whereas outside it decreases more rapidly with increasing r. For a period of maximum solar 

13 

activity. Fig. 11 of Erickson shows that N is probably about a factor of five higher in the radial 

g 

range represented by Eq. (4.2). Substituting Eq. (4.2) into (4.1), yields 


At 


m 


6.5 X 10 
f 2 d 


24 


tan_1 (-f) - tan_1 (-i)] 


sec 


(4.3) 


Here d, x^, and x are expressed in centimeters and are defined as in Fig. 1. 

and with x <0, we find 
P 


At 


m 


6.5 X 10 
f 2 d 


24 


7T 


For d « x , lx 
e 1 p 

(4.4) 


For the Arecibo Ionospheric Observatory's frequency of 430 Mcps, the lowest at which inter- 

-4 

planetary time-delay measurements are currently being made, Eq. (4.4) yields At^ ~4 x 10 sec 
for observations of Mercury near superior conjunction with d « 4R (This latter value corre- 

O 

sponds to an angular distance from the sun of 1°, the smallest at which Arecibo measurements 

-4 

can be made.) In this case, At^ would equal about 1.4 X 10 sec, as can be seen from Fig. 2, 
and would most likely be masked by the uncertainty in At . Although At^ varies inversely with 
d, whereas the corresponding dependence in At^ is logarithmic, the difference At^ — At^ is 
nowhere large enough and positive for a really reliable and accurate result to be obtained solely 
from Arecibo data. Since, for sufficiently high f, Ar^ varies as the inverse square of the radar 
frequency, this plasma effect will be reduced by a factor of almost 400 (and will therefore be 
unimportant) for measurements made at the 8350-Mcps frequency of the newly constructed, but 
not yet fully instrumented. Haystack radar facility. The Jet Propulsion Laboratory Goldstone 
radar, operated at a frequency of 2388 Mcps, could probably also make these time-delay measure¬ 
ments in a "quiet-sun" period without undue interference from the solar corona. However, during 
a maximum of solar activity, it is doubtful whether measurements at 2388 Mcps could be used 


* Note the caption on this figure: "A complication of data. ..." It unwittingly summarizes the experimental 
situation quite well. 
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alone to perform an accurate test of general relativity: The plasma would cause an increase in 

-4 -4 

time delay at d = 4R of about 0.7 x 10 sec and at d = 20R of about 0.15 x 10 sec, with an un- 

^ >]< ^ _ A 

certainty of perhaps a factor of two or three. But in the latter case, At is only 0.7 x 10 sec. 

In any event, simultaneous equivalently accurate measurements at two well-spaced frequencies 
will allow the plasma effect to be deduced and subtracted since AT m is frequency dependent and 
At is not.^ 

r 5 

Other possibly relevant effects on the delays are easily disposed of. A previous study has 

shown that the earth's and planet's atmospheres and ionospheres will not significantly affect time 

delays, even for f = 430 Mcps. The effect of the earth's gravity and motion on the laboratory 

clock is unimportant for this experiment, since the clock rate remains constant over a year to 

10 

within about one part in 10 . The gravitational effects of the earth, moon, and target planet on 

the delays are much smaller than those of the sun, but in any case the former (excepting the 
moon) could be neglected since these contributions will be almost identical in each measurement 
and consequently indistinguishable from a small decrease in the planet's radius. Similarly, any 
bias introduced into the time delay by the radar system will not affect this experiment, provided 
only that such bias is independent of the relative orientations of the earth, sun, and planet. Any 
lack of precision in the determination of c in terms of terrestrial units (such as in km/sec) is 
clearly irrelevant to our experiment since time delays only are of concern. 

In making the time-delay measurements, we must consider the radio interference introduced 
by the sun. For the Haystack facility, the antenna beamwidth is sufficiently narrow and the near 
sidelobes are of sufficiently low gain that the beam can be directed within a half-degree of the 
solar limb without the radio emanations introducing a significant increase in the over-all system 
noise temperature. For Arecibo, the closest possible approach is about 1°, partly because of 
the larger radio diameter of the sun at 430 Mcps and partly because of the higher-gain sidelobes 
of the Arecibo antenna. 

On the basis of these analyses, we can feel reasonably confident that the time-delay experi¬ 
ments discussed here will provide a meaningful test of Einstein's theory of general relativity. 
Because the magnitude of the effect decreases slowly with the distance of closest approach of the 
radar wave to the sun (logarithmic dependence on d), we have the possibility of testing not only 
a single numerical prediction, but the functional form of At^ as well. We require, of course, 

the ability to measure time delays at superior conjunction with an error of no more than about 

• \ 4 

lOpsec. The upgraded Haystack facility is expected to provide such a capability. Repeated 

measurements extended over a period of several years should then enable the maximum effect 

to be determined to about 1 percent and the effect at greater distances of closest approach to 

correspondingly lower accuracies (see Fig. 2). We need not worry about constant biases in the 

radar system limiting our use of statistics to improve accuracy since, as shown above, constant 

biases serve only to introduce an error in the estimate of the planet radius. 

In principle, this relativistic effect on time delay could also be observed either by using space 

probes in orbit about the sun, or by placing transponders on a planet. At present, each of these 

approaches would be more difficult and more costly to implement than the one discussed here. 

* Since for r > 20R $ , N(r) decreases more rapidly with increasing r than is indicated in Eq. (4. 2), we might find 
AT m for d = 20R S to be somewhat less than estimated. 

t It may be possible to estimate the plasma effect on delay adequately by simultaneously observing the deflection 
of low-frequency (<100 Mcps) signals from suitable radio stars near the Iine-of-sight to the target planet. How¬ 
ever, a variety of plasma conditions having different effects on delay can lead to the same over-all deflection of 
the radio star emanations. 
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V. ADDITIONAL SCIENTIFIC "FALLOUT" 


The time-delay measurements proposed to test the predictions of general relativity con¬ 
cerning the effect of gravitational fields on light rays also have other scientific implications. 

The orbital accuracies achievable will represent not only a very significant improvement over 
present knowledge, but will enable a somewhat independent and more accurate determination to 
be made of the precession of Mercury’s perihelion position. With the measurements distributed 
fairly uniformly about the orbits of both earth and Mercury, the standard deviation a(cj) of the 
error in perihelion precession is given by 


a(d>) 


2 X 10 
(Nt 3 ) 1/2 


— sec of arc/100 yr 

ae t . r j 


min 


(5.1) 


where N is the number of measurements made per year, t is the total time span of the observa¬ 
tions in years, g(t) is the standard deviation of the single-measurement error, r ^ is the mini¬ 
mum range to Mercury, and a (in astronomical units) and e are the semimajor axis and eccen¬ 
tricity, respectively, of Mercury’s orbit. (The numerical coefficient merely represents the 
conversion from radians per revolution to seconds of arc per century.) With N = 50, t = 3, 

c t(t) = 10 usee, t . = 350 sec, a = 0.38, and e = 0.2, we find 

' r mm 

u(cj) « 0.2 sec of arc/100 yr , (5.2) 

which is considerably smaller than the error currently attributed to the conventional determina¬ 
tion of the centennial value of the precession. 

In order to determine to this accuracy the anomalous (non-Newtonian) contribution to the 
precession, we must of course know the Newtonian contribution to a comparable degree of ac¬ 
curacy. Two major impediments to this determination are the uncertainty in our knowledge of 
the sun’s quadrupole moment and of Venus's mass. The former problem has been discussed 
elsewhere.^ The latter can probably be solved by using earth-Venus time-delay observations 
from which the perturbations of Venus on the orbit of earth should enable Venus’s mass to be 
determined accurately enough. Knowledge of the mass of Mercury may also be improved sig¬ 
nificantly in virtue of its effect on the orbits of earth and Venus. 

The determination of the radii of both Venus and Mercury will also be refined substantially, 
in fact by about two orders of magnitude. If the optical diameter of Venus were determinable 
with an accuracy of about ±5 km, then the height above the surface of the Venusian cloud layer 
could be deduced with reasonable accuracy; such a result would be of great importance for the 
study of planetary atmospheres.t The combination of accurate radii and mass determinations 
will, of course, yield correspondingly accurate densities for Mercury and Venus. These values 
are of interest to studies of planetary formation. 


* Although it is claimed by Jet Propulsion Laboratory analysts that their tracking data on Mariner li yielded five- 
place accuracy in the determination of the mass of Venus, there is the possibility that systematic errors may result 
in this accuracy being more apparent than real. 

t There is also the possibility that the X-band radar waves are reflected from a cloud layer and the S-band and 
lower frequency waves from the surface. In such a case, simultaneous delay measurements at X-band and at a 
lower frequency might enable the height of the reflecting layer to be determined. 
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Highly accurate determinations of the mean anomalies of the inner planets will allow a rea¬ 
sonably strict limit to be placed on the possible time dependence of the gravitational constant G. 

17 

[The theoretical conjecture (first made by Dirac ) that G may decrease with time is based 
partly on the apparent expansion of the universe.] If the mean anomalies could be determined 
each year with an error no greater than 0.02 seconds of heliocentric arc, which should be achiev¬ 
able if individual time-delay measurements are made with errors of lOpsec or less, then a var- 

10 

iation in G of about five parts in 10 per year would be discernible after several years of ob- 

>Jc 

servations. Of course, one must also have clocks with the requisite long- and short-term 

stabilities. Since the theoretical estimates of the rate of change of G vary down to a part in 
11 18 

10 per year (see, e.g., Dicke ), this test will probably not be crucial unless the measure¬ 
ments are either improved in accuracy or continued over many years. 

As stated in Sec. IV, the solar corona affects time-delay measurements by an amount in¬ 
versely proportional to the square of the radar frequency. In particular, accurate delay measure¬ 
ments made at about 400 Mcps, for example at Arecibo, would afford an excellent opportunity 
to study the solar corona. The integrated electron densities could be determined as a function of 
d and from these results the average radial electron density could be deduced. Short-term 
fluctuations could be studied by making frequent measurements and long-term trends observed 
by extending the study over the length of the solar cycle. The solar-corona plasma will also 
influence the Doppler shift of the radar waves via an effect analogous to the one discussed at the 
end of Sec. III. As in that discussion, we conclude that the time-delay differential holds more 
promise at present as a tool for investigating the corona. However, Doppler observations may 
enable the detection of dense plasma "wedges" moving perpendicular to the path of the radar 
wave. Bending caused by the plasma would be significant only close to the solar limb where 
measurements cannot be made. Hence this effect too seems to be of less interest than the time- 
delay measurements. 

Faraday rotation effects, caused by the different propagation speeds of the ordinary and 
extraordinary modes, would also enable certain characteristics of the solar magnetic field to be 
investigated. The detailed deductions are somewhat complicated, however, because sign re¬ 
versals in the magnetic field tend to cancel in their effect on Faraday rotation. In summary, 
such radar experiments could probably provide all the solar corona information currently ex¬ 
pected from the forthcoming Sun-Blazer probes. 


* Note that the effect of a change in G on the mean anomaly is proportional to the square of the elapsed time, 
whereas the corresponding effect on the Astronomical Unit is merely proportional to the time. (In this discussion, 
we are, of course, referring to nongravitational (e.g., atomic) measures of time.) 
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